Transition from a two-dimensional superfiuid to a one-dimensional Mott insulator 
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A two-dimensional system of atoms in an anisotropic optical lattice is studied theoretically. If 
the system is finite in one direction, it is shown to exhibit a transition between a two-dimensional 
superfiuid and a one-dimensional Mott insulating chain of superfiuid tubes. Monte Carlo (MC) 
simulations are consistent with the expectation that the phase transition is of Kosterlitz-Thouless 
(KT) type. The effect of the transition on experimental time-of-flight images is discussed. 



Cold atomic gases in optical lattices provide a means to 
study many-body quantum phenomena that offers both 
versatility, precision, and control. These techniques al- 
lowed for the spectacular realization of the Mott tran- 
sition by Greiner et al. ,l|, which has inspired numer- 
ous proposals for creating and detecting exotic quantum 
phases in such systems. 

By manipulating the alignment and the polarizations 
of the laser beams that build up the optical lattice, the 
geometry of the lattice can be changed. The probability 
for tunneling can be adjusted by changing the separa- 
tion between potential wells, or by simply varying the 
irradiancc. In conjunction, this brings about the possi- 
bility of designing optical lattices where the probability 
for tunneling differs significantly between different direc- 
tions, and where this difference in tunneling rates can be 
a control parameter. By allowing the atoms to tunnel 
only along one direction, one essentially has created a 
2D lattice of ID quantum gases. By allowing tunneling 
in two directions, an array of 2D systems is created 0. 

The tunability allows to explore the crossover between 
different dimensionalities. Thus, with a properly chosen 
geometry, it should be possible to start off with high laser 
irradiance and thus with deep potential wells. This sys- 
tem will be a Mott insulator, where the number of par- 
ticles per site is fixed to an integer and there is no phase 
coherence. By lowering the irradiance in one direction, 
increasing the tunneling probability, one should reach a 
point where tunneling becomes probable along just one 
direction. This would give us isolated ID tubes, known 
as Luttinger liquids [3(. Decreasing the irradiance in a 
second direction, we expect to reach a point where there 
is a crossover from this to a system of 2D superfluids, 
and eventually one global 3D superfiuid 0, [B|, Q • In this 
Letter, this dimensional crossover is addressed by quan- 
tum MC simulations of an anisotropic bosonic Hubbard 
model in two dimensions. We show how the model can be 
readily realized in experiment, and predict the outcome 
of absorption imaging on both sides of the transition. 

System. - Consider a gas of bosonic atoms at zero 
temperature in a 2D, anisotropic optical lattice. The gas 



is described by the Bose-Hubbard Hamiltonian Q 



U 



— ^2ni(m - 1) - mX!" 1 ' 



(1) 



where i = (i x ,i y ); bi, b] are bosonic annihilation and cre- 
ation operators, = b\bi, U the interaction strength, /i 
the chemical potential, and t x , t y the tunneling matrix 
elements in the Cartesian directions. The crossover from 
2D to ID behavior considered in this paper can be re- 
alized in a simple cubic optical lattice with a potential 
of the form V(r) = J2 a=x y z^ a cos2 ^ tnc poten- 
tial barrier in the z direction V z is strong enough, there 
can be no tunneling in the z direction, and the sample 
can be considered as a stack of independent 2D systems 
. V x should be chosen relatively weak in order to allow 
for superfiuid flow in the x direction, and V y is scanned 
over the physically interesting range. As will be shown, 
if the lattice has an extent of eight sites in the x direc- 
tion there will occur a phase transition for t x /U = 0.3 
and t y /U » 6 x 10~ 3 . As an example, for a gas of 87 Rb 
atoms in an optical lattice with wavelength A = 850 nm, 
these tunneling matrix elements are obtained by choos- 
ing potential strengths V x = 2.5i? rcc , V y — 3QE ICC , and 
V z = 80-E rcc , where E lec is the single-photon recoil energy 
for a wavelength of 780nm 

Phase diagram. - The Hubbard model exhibits a zero- 
temperature quantum phase transition from a superfiuid 
to a Mott insulating state when the ratio of the tunnel- 
ing matrix elements and the coupling strength, t x /U and 



ty/U, are decreased below a critical value [8j. In one 
dimension (i.e., either t x = or t y = 0), it is known 
that the phase transition occurs at t/U = 0.3 for the 
system with an average of one particle per site [9(. For 
the anisotropic 2D model considered here, a mean-field 
argument implies that t should be replaced by the sum 
t x +ty in the ID expression, so the transition occurs when 
(t x + t y )/U = 0.3. For finite systems the phase transition 
occurs at a smaller value of t/U. 

Consider the anisotropic case, where one of the tun- 
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FIG. 1: Phase diagram for the finite anisotropic Hubbard 
model at number density n = 1. The 2D superfluid and Mott 
insulating phases are denoted by "2D SF" and "2D MI", re- 
spectively. The size of the ID Mott insulating phase "ID MI" 
has been exaggerated in the figure. The vertical double arrow 
indicates the interval along which the MC calculations in this 
study have been performed. 
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FIG. 2: Velocity distribution along the weakly coupled lattice 
direction. The curves corresponds to a ratio of the tunnel- 
ing matrix element to the interaction strength, t y /U = 0.012 
(dashed line), 0.0025 (full), and 1.5 x 10~ 3 (dotted). The inset 
shows the full 2D profile of the velocity distribution for the 
case ty/U — 1.5 x 10 -3 . The system size is 12x12 sites. 



neling matrix elements, say, t x , is larger than the critical 
ID value and the other one, t y , is varied. In an infi- 
nite 2D system, it is known that the system remains a 
2D superfluid (2D SF) right down to the point where 
t y vanishes [Toj . In other words, when the coupling in 
one direction vanishes, the system can be seen as an ar- 
ray of uncoupled Luttinger liquids. However, any finite 
tunneling, no matter how small, will induce phase co- 
herence between the Luttinger liquid tubes, so that an 
anisotropic phase, with superflow only in one direction, 
does not exist. However, the situation changes when the 
system is not infinite. If each Luttinger liquid has a fi- 
nite length, there is an energy barrier against tunneling 
into and out of it, and as a result there exists, for small 
t y , an anisotropic state with coherence only along the x 
direction @, d, @]. Hence, when the weaker tunneling 
matrix element t y is small enough, the state of the sys- 
tem is a ID Mott insulator (ID MI), where each site is a 
Luttinger liquid tube extending along the x direction. If 
the tunneling t y is increased above a critical value t c , the 
tubes become phase correlated and the system undergoes 
a quantum phase transition to a superfluid state. In a 
geometry where the system is infinite in the weakly cou- 
pled direction, the state is an ID SF, i. e., a superfluid 
chain of Luttinger liquids, just above the phase transi- 
tion, but is expected to undergo a crossover into a 2D SF 
as the coherence length becomes smaller than the tube 
length L x , and the picture of separate Luttinger liquid 
tubes can no longer be sustained. The crossover between 
ID SF and 2D SF presumably happens over a very short 
interval, why we in the following shall speak of the ID 
MI - 2D SF transition. In an actual experiment where 
the system is finite, there is a crossover between ID MI 
and 2D SF without an intermediate ID SF state. In the 



ID MI phase, each tube along the x direction behaves 
as an independent Luttinger liquid, characterized by an 
(inverse) exponent K [!, 0, [|| . This so-called Luttinger 
parameter determines the algebraic decay of the single- 
particle correlation function according to 

r x (j)^(bl tiy b ix+j , iy )~ r WK*). ( 2) 

This form for the asymptotic behavior holds only in the 
ID case, i. e., for decoupled Luttinger liquid tubes. T y is 
defined analogously. The critical tunneling was in Refs. 
0, S S| found to depend on the parameters as 

t c otL~ 2+ ^. (3) 

Since it is known that K x > 2 here Q, the exponent in 
Eq. is expected to be slightly less than -1.75. The 
constant of proportionality in Eq. ([3]) depends on the 
model used to integrate out the dynamics in the tubes 
as well as on the dimensionality. These phase boundaries 
have been indicated at arbitrarily chosen positions in the 
phase diagram in Fig. [TJ 

Experimental signal. - Detection of the transition can 
be achieved by applying an optical lattice potential, V(r), 
to a Bose-Einstein condensate in a magnetic trap. Ramp- 
ing up V(r) adiabatically ensures that the atoms are not 
heated, and that they are loaded to the lowest energy 
band of the lattice [l[. For high lattice barriers, the cor- 
relations between wells are lost after some hold time, re- 
sulting in a Mott insulating state, whereas in the case of 
low barrier heights the sample stays coherent. The state 
of the system is detected by standard time of flight ab- 
sorption imaging. When the sample is released by ramp- 
ing down the lattice and magnetic trap non-adiabatically, 



it will expand freely. The expansion reveals the sample's 
velocity distribution n(k), given by the Fourier trans- 
form of the trapped sample density distribution in real 
space n(r) n(k) oc |u>(k)| 2 5(k), where the first part 
is the Fourier transform of the Wannier function for the 
first Bloch band w(r), and the second part, the structure 
factor S(k), is the Fourier transform of the correlation 
function T(r) defined in Eq. |(5J). The structure factor 
has the periodicity Air/\. 

Knowing the structure factor from MC simulations (see 
below), and calculating the Wannier functions for a cu- 
bic lattice, a theoretical prediction for the actual exper- 
imentally observable time-of-flight profiles can thus be 
obtained. The result is shown in Fig. [2l The profile 
of the cloud in the direction of weak tunneling is deter- 
mined by the state of the system. When distinct peaks 
are visible at k = 0, ±4ir/\, a 2D SF is detected. This is 
connected with a slow decay of the real-space correlation 
function T(r). If the lattice depth is increased, the sam- 
ple passes into the ID MI and the interference fringes are 
blurred as T(r) becomes exponentially decaying. The sig- 
nal would be even more pronounced for a larger system. 
However, the peaks will be visible far into the Mott insu- 
lator regime, an effect caused by maintained short-range 
coherence [13(. 

Method. - To solve for the ground state of the Bose- 
Hubbard model, Eq. {l}, Quantum MC calculations with 



the stochastic series expansion algorithm are used [11 
12j . We simulate a system of size L x x L y lattice sites with 
periodic boundary conditions to minimize edge effects. 
The inverse temperature is (3 = (fcsT) -1 = 1000C/ -1 ; 
this is large enough to ensure that the system is in its 
ground state. The results presented in the article are 
calculated with a fixed t x = 0.3U with t y /U being varied 
across the anticipated phase transition as indicated by 
the double arrow in Fig. [TJ Trial calculations not shown 
here, using larger values of t x /U , have yielded similar 
results. The calculations have been performed at the 
fixed average number of particles per site n = 1. 

Transition. - To locate the 2D SF<-» ID MI transition, 
three different observable quantities are calculated, each 
of which can be considered a measure of superfluidity. 
The superfluid density, p s , is readily obtained from the 
simulation data as p s = (Wy)/f3, where W 2 is the so- 
called square winding number in the y direction, defined 
as the squared net number of times a particle crosses 
the periodic boundary in the calculations [3]. Second, 
the off-diagonal correlation function along the y direc- 
tion, Ty(J), is, as we have seen, experimentally accessible 
through the momentum distribution, which is closely re- 
lated to its Fourier transform. We now study the corre- 
lation function in real space evaluated at half the system 
size, T y {y = L y /2) which we expect to behave differently 
above and below the transition. Above, T y (y) decays 
algebraically, below exponentially. A third observable 
is the compressibility, defined as the mean of variances, 
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FIG. 3: From top to bottom, the density fluctuation, the cor- 
relation, and the superfluid density are shown as a function of 
the size-independent parameter tyL^' 75 /U, for different sys- 
tem sizes L x x L y measured in lattice sites, as indicated in the 
legend. The scale is logarithmic on both axes. The correlation 
function and the compressibility are multiplied by arbitrarily 
chosen numerical factors to separate the three groups of lines. 
The superfluid density is divided by the tunneling matrix el- 
ement for better visual clarity. 



AJV 2 = Yli var (Si x n i*,i y ) /Ly It is expected that 
AN 2 changes its behavior at the transition: due to the 
finite length of the tubes, the fluctuations in the number 
of particles per tube will be suppressed once the system 
becomes insulating in the weakly coupled y direction. 

In Fig. [31 the three quantities, p s , T y (L y /2), and AA 2 , 
are shown for lattices 6 < L x = L y < 14. The observ- 
ables are plotted as functions of i y Lj. 75 /U . The exponent 
of L x is given by Eq. ([3]) and the value of K x is deter- 
mined by a finite-size scaling of the superfluid density 
as explained below. For a system with t x /U — 0.4, not 
shown here, the best scaling is obtained for t y L x 79 /U. 
The superfluid density in the y direction is in the super- 
fluid phase seen to be approximately proportional to the 
coupling constant t y /U, hence the superfluid density is 
divided by t y /U. It is clear from Fig. [3] that there is a 
change in behavior around t y /U » 0.3L~ 175 . The three 
quantities in Fig. [3J are behaving as in a ID MI-SF tran- 
sition. The true ID system can be obtained by putting 
L x = 1 and letting L y — > oo, and we then recover the 
known value for the true ID transition t y /U — 0.3 0]. 
(The same curves are obtained for L y > L x , but these 
data are not shown here not to overload the figure). The 
finite size effects above the transition point are taken care 
of by scaling L x with the power 1.75. If the curves in- 
stead were plotted as a function of t y /U, they would not 
tend to coincide. 

Since in the critical region, the system can be described 
as a ID chain of Luttinger liquid tubes, the phase transi- 
tion is a ID quantum phase transition. Such a transition 
can be mapped onto a classical phase transition in 2D 
[l5| , and if the number of particles is held fixed as in the 
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FIG. 4: a) Superfluid density in the dual XY model, Eq. (TJJ, 
scaled according to Weber-Minnhagen scaling, b) Exponent 
of the algebraic decay of the correlation Ty, as a function of 
the dimensionless coupling parameter tyLl.' 75 /U . The hori- 
zontal line shows the critical value expected for a ID Mott 
transition at fixed density, 1/(2K) — 0.25. Different curves 
correspond to different system sizes; the same symbols are 
used in both parts of the figure. 



present study, it is known that the transition considered 
here is a KT transition Q. 

Characteristic for such a transition is a discontinuous 
jump of the superfluid density at the critical point [14j . In 
Fig.|3]a) a finite size scaling of the data for the superfluid 
density is illustrated. In order to do the scaling, a series 
of mappings has to be performed: First, according to 
Ref. [4[ , the original Hubbard model can be rewritten as 
a Josephson junction chain model, where each Luttinger 
liquid tube is treated as a site, with a charging energy 
Ec oc L" 1 and hopping energy Ej oc t y L x 1 / 2K:c _ Next, 
this model is in turn mapped onto a 2D XY-model if 
one identifies f3 X Y = \fEj/E c , L x . X y = (5\/EjEc, and 
L y ,XY = L y |15j |. where the quantities referring to the 
dual 2D XY-model carry the subscript XY. As a result, 
the superfluid density in the 2D XY-model is [3] 



Ps,XY 



(X (W^Ly/PtyL 1 - 1 / 2 ^ 



(4) 



The quantum phase transition occurring at a critical tun- 
neling in the underlying Hubbard model corresponds to 
the well-known finite-temperature KT transition in the 
2D XY-model. According to Weber-Minnhagen scaling 
[3) Ps,xy should be proportional to 1 + 1/(2 In i y + C) 
at the critical point, where C is a fitting constant to be 
determined. The quantity p s ,XY divided by this func- 
tion should therefore assume the same value for all sys- 
tem sizes at the critical point, if C and K x are chosen 
correctly. The best scaling is obtained with C = — 1 and 
K x — 2 and is shown in Fig. 0] a). The scaling behav- 
ior supports the conclusion that the transition is a KT 
transition and the critical point is t c /U « 0.27L" 175 . 

An independent method for calculating the location of 
the transition is provided by fitting the correlation func- 
tion in the y-direction, T yi to an algebraic decay func- 
tion according to Eq. @. The phase transition for a ID 



Hubbard model at a fixed number of atoms per site is 
found by determining the point where the exponent in 
the power law, l/(2K y ), becomes equal to 0.25 Q. 

In Fig. 2] b) the exponent obtained from an algebraic 
fit is displayed for different system sizes. There is a clear 
finite-size effect, the larger L y the larger the critical value. 
The critical point at which l/(2K y ) = 0.25 for the infi- 
nite system is determined to be t c /U « 0.3L^ 175 , in 
agreement with the scaling in Fig. [5] a) and the point of 
loss of superfluidity observed in Fig. [3] 

Summary - We have reported theoretical evidence for 
a transition in a 2D Hubbard model between a 2D super- 
fluid and a ID Mott insulating state consisting of isolated 
ID tubes. The location of the transition is consistent 
with predictions based upon a random phase approxi- 
mation made in Rcfs. 

@, H, II- By scaling it is shown 
that the system undergoes a KT transition. The transi- 
tion can be induced using atoms in an optical lattice and 
observed by time-of-flight imaging. 
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